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i Abstract 



Wc prove that the pleated hyperbolic paraboloid, a familiar origami model known since 1927, 
in fact cannot be folded with the standard crease pattern in the standard mathematical model 
' of zero-thickness paper. In contrast, we show that the model can be folded with additional 

creases, suggesting that real paper "folds" into this model via small such creases. We conjecture 
that the circular version of this model, consisting simply of concentric circular creases, also folds 
without extra creases. 

' At the heart of our results is a new structural theorem characterizing uncreased intrinsically 

t/3 , flat surfaces — the portions of paper between the creases. Differential geometry has much to 

say about the local behavior of such surfaces when they are sufficiently smooth, e.g., that they 
are torsal ruled. But this classic result is simply false in the context of the whole surface. Our 
structural characterization tells the whole story, and even applies to surfaces with discontinuities 
^ , in the second derivative. We use our theorem to prove fundamental properties about how paper 

folds, for example, that straight creases on the piece of paper must remain piecewise-straight 
. (polygonal) by folding. 

^ '. 1 Introduction 

o\ 

. A fascinating family of pleated origami models use extremely simple crease patterns — repeated 

concentric shapes, alternating mountain and valley — yet automatically fold into interesting 3D 
k> I shapes. The most well-known is the pleated hyperbolic paraboloid, shown in Figure [H where the 

j_j ■ crease pattern is concentric squares and their diagonals. As the name suggests, it has long been 



conjectured, but never formally established, that this model approximates a hyperbolic paraboloid. 
More impressive (but somewhat harder to fold) is the circular pleat, shown in Figure [2l where the 
crease pattern is simply concentric circles, with a circular hole cut out of the center. Both of these 
models date back to the Bauhaus, from a preliminary course in paper study taught by Josef Albers 
in 1927-1928 |Win69L p. 434], and taught again later at Black Mountain College in 1937-1938 
|Adl04l pp. 33, 73]; see |DD08j . These models owe their popularity today to origamist Thoki Yenn, 
who started distributing the model sometime before 1989. Examples of their use and extension for 
algorithmic sculpture include [DDL99l[K"DD08j . 
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(a) Standard mountain-valley pattern. (b) Photograph of physical model. [Jenna Fizel] 

Figure 1: Pleated hyperbolic paraboloid. 




(a) Mountain-valley pattern. (b) Photograph of physical model. [Jenna Fizel] 

Figure 2: Circular pleat. 



The magic of these models is that most of the actual folding happens by the physics of paper 
itself; the origamist simply puts all the creases in and lets go. Paper is normally elastic: try 
wrapping a paper sheet around a cylinder, and then letting go — it returns to its original state. But 
creases plastically deform the paper beyond its yield point, effectively resetting the elastic memory 
of paper to a nonzero angle. Try creasing a paper sheet and then letting go — it stays folded at the 
crease. The harder you press the crease, the larger the desired fold angle. What happens in the 
pleated origami models is that the paper tries to stay flat in the uncreased portions, while trying 
to stay folded at the creases, and physics computes a configuration that balances these forces in 
equilibrium (with locally minimum free energy). 
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But some mathematical origamists have wondered over the years |Wer05j : do these models 
actually exist? Is it really possible to fold a piece of paper along exactly the creases in the crease 
pattern of Figures [1] and [2p The first two authors have always suspected that both models existed, 
or at least that one existed if and only if the other did. But we were wrong. 



Our results. We prove that the hyperbolic-paraboloid crease pattern of Figure 1(a) does not 
fold using exactly the given creases, even with a hole cut out of the center. 

In proving the impossibility of folding the pleated hyperbolic paraboloid, we develop a structural 
characterization of how uncreased paper can fold (hence the title of this paper). Surprisingly, such 
a characterization has not been obtained before. An intuitive understanding (often misquoted) is 
that paper folds like a ruled surface, but that claim is only true locally (infinitesimally) about every 
point. When the paper is not smooth or has zero principal curvature at some points, the truth gets 
even subtler. We correct both of these misunderstandings by handling nonsmooth (but uncreased) 
surfaces, and by stating a local structure theorem flexible enough to handle zero curvatures and all 
other edge cases of uncreased surfaces. 

In contrast, we conjecture that the circular-pleat crease pattern of Figure 2(a) folds using exactly 
the given creases, when there is a hole cut out of the center. A proof of this would be the first 
proof of existence of a curved-crease origami model (with more than one crease) of which we are 
aware. Existing work characterizes the local folding behavior in a narrow strip around a curved 
crease, and the challenge is to extend this local study to a globally consistent folding of the entire 
crease pattern. 

Another natural remaining question is what actually happens to a real pleated hyperbolic 
paraboloid like Figure 1(b), One conjecture is that the paper uses extra creases (discontinuities in 
the first derivative), possibly many very small ones. We prove that, indeed, simply triangulating 
the crease pattern, and replacing the four central triangles with just two triangles, results in a 
foldable crease pattern. Our proof of this result is quite different in character, in that it is purely 
computational instead of analytical. We use interval arithmetic to establish with certainty that 
the exact object exists for many parameter values, and its coordinates could even be expressed by 
radical expressions in principle, but we are able only to compute arbitrarily close approximations. 



2 Structure of Uncreased Flat Surfaces 

Our impossibility result rests on an understanding of how it is possible to fold the faces of the 
crease pattern, which by definition are regions folded without creases. The geometric crux of the 
proof therefore relies on a study of uncreased intrinsically flat (paper) surfaces. This section gives 
a detailed analysis of such surfaces. Our analysis allows discontinuities all the way down to the 
second derivative (but not the first derivative — those are creases), provided those discontinuities 
are somewhat tame. 

We begin with some definitions, in particular to nail down the notion of creases. 

Definition 1 For us a surface is a compact 2-manifold embedded in M^. The surface is if 
the manifold and its embedding are . The surface is piecewise-C^ if it can be decomposed as a 
complex of regions joined by vertices and edges. 

Definition 2 A good surface is a piecewise-C^ surface. A good surface S therefore decomposes 
into a union of C'^ surfaces Si, called pieces, which share edges 7j, called semicreases, whose 
endpoints are semivertices. Isolated points of discontinuities are also semivertices. If S is itself 
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everywhere on a semicrease, we call it a proper semicrease; otherwise it is a crease. Similarly a 
semivertex v is a vertex if S is not at v. Accordingly an uncreased surface is a good surface 
(with no creases or vertices), and a creased surface is a good surface not everywhere (with at 
least one crease or vertex). 

Definition 3 A surface is (intrinsically) flat if every point p has a neighborhood isometric to a 
region in the planea 

In order to understand the uncreased flat surfaces that are our chief concern, we study the 
flat surfaces that make them up. On a surface, the well-known principal curvatures ki > K2 
are defined for each interior point as the maximum and minimum (signed) curvatures for geodesies 
through the point. A consequence of Gauss's celebrated Theorema Egregium |Gau02j is that, on a 

flat surface, the Gaussian curvature k\K2 must be everywhere zero. Thus every interior point 
of a flat surface is either parabolic with /c2 7^ ^1 = or planar with k2 = ki = 0. 

Each interior point p on a flat surface therefore either 

(a) is planar, with a planar neighborhood; 

(b) is planar and the limit of parabolic points; or 

(c) is parabolic, and has a parabolic neighborhood by continuity, 
and an interior point on an uncreased flat surface may additionally 

(d) lie on the interior of a semicrease; or 

(e) (a priori) be a semivertex. 

For points of type (a), it follows by integration that the neighborhood has a constant tangent 
plane and indeed lies in this plane. Types (b) and (c) are a bit more work to classify, but the 
necessary facts are set forth by Spivak |Spi79[ vol. 3, chap. 5, pp. 349-362] and recounted below. (In 
Spivak's treatment the regularity condition is left unspecified, but the proofs go through assuming 
only C^.) We address type (d) farther below. From our results it will become clear that the 
hypothetical type (e) does not occur in uncreased flat surfaces. 



Proposition 1 |Spi79 Proposition III. 5. 4 et seq.] For every point p of type (c) on a surface M , 



a neighborhood U C M of p may be parametrized as 

f{s,t) = c{s)+t-6{s) 

where c and 6 are functions; c(0) = p; \5{s)\ = 1; c'(s), 5{s), and 5'{s) are coplanar for all s; 
and every point of U is parabolic. 

Write interior(M) for the interior of a surface M. 

Proposition 2 |Spi79[ Corollaries III. 5. 6-7] For every point p of type (b) or (c) on a surface M, 
there is a unique line Lp passing through p such that the intersection LpCi M is open in Lp at p. 
The component Cp containing p of the intersection Lp n interior (M) is an open segment, and every 
point in Cp is also of type (b) or (c) respectively. 



^Henceforth we use the term "flat" for this intrinsic notion of the surface metric, and the term "planar" for the 
extrinsic notion of (at least locally) lying in a 3D plane. 
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Following the literature on flat surfaces, we speak of a segment like the Cp of Proposition [2] as 
a rule segment. The ruling of a surface is the family of rule segments of all surface points, whose 
union equal the surface. A ruling is tarsal if all points along each rule segment have a common 
tangent plane. 

To characterize points of type (d), lying on semicreases, we require the following two proposi- 
tions. 

Proposition 3 Consider a point q of type (d) on a surface M . Then q is not the endpoint of the 
rule segment Cp for any point p £ Ad of type (h) or (c). 

Proof: It suffices to show the conclusion for p of type (c), because a rule segment of type (b) is a 
limit of rule segments of type (c). Let 7 be the interior of the semicrease on which q lies. 

Because M is C^, it has a tangent plane Mq at each q, which is common to the two pieces 
bounded by 7. Parametrize 7 by arclength with 7(0) = q, and write n(s) for the unit normal to 
the tangent plane M^^^^y Parametrize the two pieces as torsal ruled surfaces by the common curve 
ci(s) = C2(s) = 7(5) and lines 5i{s) and 52{s). Then, because each piece is torsal, n(s) _L 5i(s) and 
n{s) _L 52{s). But both 5i{s) and 52{s) lie in the tangent plane at s, perpendicular to n(s), and so 
too does n(s) because n(s) is always a unit vector. Therefore, for each s, either h{s) = or 6i{s) 
and (52 (s) are collinear. 

Let A be the subset of 7 on which n(s) = 0, and B the subset on which 6i{s) and ^2(5) are 
collinear. Then we have shown that AuB = 7. By continuity, both A and B are closed. Therefore 
any point of 7 which does not belong to A is in the interior of i?, and any point not in the interior 
of A is in the closure of the interior of B. 

If an open interval / along 7 is contained in A so that h{s) = 0, then a neighborhood in M of I 
is planar by integration because each rule segment has a single common tangent plane in a torsal 
ruled surface. On the other hand if / is contained in B so that 5i{s) and 62{s) are collinear, then a 
neighborhood is a single ruled surface. In either case, the rule segments from one surface that 
meet / continue into the other surface. That is, each rule segment meeting a point in the interior 
oi A or B continues into the other surface. 

Now we conclude. By continuity, each rule segment meeting a point in the closure of the interior 
of A or continues into the other surface; but these two closures cover 7. So no rule segment ends 
on 7, including at g. □ 



Proposition 4 For every point p of type (d) on a surface M , there is a unique line Lp passing 
through p such that the intersection LpD M is open in Lp at p. The component Cp containing p 
of the intersection Lp D interior (M) is an open segment, the limit of rule segments through points 
neighboring p, and every point of Cp is also of type (d). 

Proof: Let Br{p) be a radius-r disk in M centered at p, small enough that no point of the disk 
is a semivertex. By Proposition [3l the rule segment through any point q of type (b) or (c) in the 
half-size disk Br/2{p) cannot end in Br{p), so it must be of length at least r/2 in each direction. 
Further, p must be a limit of such points, or else a neighborhood oi p would be planar. 

By a simple compactness argument, provided in |Spi79| for the type-(b) case of Proposition[2l Cp 
is the limit of (a subsequence of) rule segments Cq through points q of type (b) and (c) approaching 
p and is an open segment. Because each Cq has a single tangent plane, the discontinuity in second 
derivatives found at p is shared along Cp. □ 

Two corollaries follow immediately from Proposition [H 
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Corollary 5 Every (proper) semicrease in an uncreased flat surface is a line segment, and its 
endpoints are boundary points of the surface. 

Corollary 6 An uncreased flat surface has no interior semivertices; every interior point is in the 
interior of a C"^ piece or a semicrease. 

Another corollary summarizes much of Propositions [2] and S] combined. 

Corollary 7 Every interior point p of an uncreased flat surface M not belonging to a planar 
neighborhood belongs to a unique rule segment Cp. The rule segment's endpoints are on the boundary 
of M, and every interior point of Cp is of the same type (h), (c), or (d). 

Finally, we unify the treatment of all types of points in the following structure theorem for 
uncreased flat surfaces. The theorem is similar to Proposition [H which concerns only points of 
type (c). 

Theorem 8 Every interior point of an uncreased flat surface has a neighborhood that is a ruled 
surface. In each rule segment, every interior point is of the same type (a), (b), (c), or (d). The 
ruled surface may be parametrized as 

f{s,t) = c{s)+t-6{s), 
where c is , 5 is C", and S is whenever c{s) is of type (a), (b), or (c). 

Proof: Let p be a point on an uncreased flat surface M. lip is of type (a), then we may parametrize 
its planar neighborhood as a ruled surface almost arbitrarily. Otherwise, p is of type (b), (c), or 
(d) and has a unique rule segment Cp. 

Embed a neighborhood U C M of Cj, isometrically in the plane, by a map (p : U — > M^. Let 
7 be a line segment in the plane perpendicularly bisecting (j)(Cp), parametrized by arclength with 
7(0) = (pip). Every point (/)~^{'y{s)) of type (b), (c), or (d) has a unique rule segment C^-i(^(^)); 
for such s, let e{s) be the unit vector pointing along i?!'(C(^-i('y(s))), picking a consistent orientation. 

Now the remaining s are those for which (p~^{'y{.s)) is of type (a). These s form an open 
subset, so that for each such s there is a previous and next s not of type (a). For each such s, we 
can determine an e(s) by interpolating angles linearly between the e(s) defined for the previous 
and next s not of type (a). The resulting function e{s) is continuous and identifies a segment 
through every point in 7, giving a parametrization of a neighborhood of 7 as a ruled surface by 
g{s,t) = j{s)+t-e{s). 

Finally, write /(s,t) = cl)~^{g{s,t)) to complete the construction. □ 

3 How Polygonal Faces Fold 

If all edges of the crease pattern are straight, every face of the crease pattern is a polygon. We 
first show that, if the edges of such a polygon remain straight (or even piecewise straight) in space, 
then the faces must remain planar. 

Theorem 9 If the boundary of an uncreased flat surface M is piecewise linear in space, then M 
lies in a plane. 
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Proof: Let p be a parabolic point in the interior of M, a point of type (c). We will show a 
contradiction. It then follows that every point of M is of type (a), (b), or (d), so planar points of 
type (a) are dense and by integration M lies in a plane. 

Because p is parabolic, it has by Proposition [T] a neighborhood consisting of parabolic points 
which is a ruled surface. By Corollary [71 the rule segment through each point in this neighborhood 
can be extended to the boundary of M. Let U be the neighborhood so extended. 

Now the boundary of U consists of a first rule segment ab, a last rule segment cd, and arcs bd 
and ac of which at least one must be nontrivial, say bd. Because we extended U to the boundary 
of M and the boundary of M is piecewise linear, bd consists of a chain of segments. Let b'd' be one 
of these segments. 

Let q be any point interior to the segment b'd', and consider the normal vector n{q) to M 
at q. The normal is perpendicular to b'd' and to the rule segment Cq meeting q. Because U is 
torsal, its derivative n'(q) along b'd' is also perpendicular to Cq, and because the normal is always 
perpendicular to b'd' the derivative is perpendicular to b'd'. But this forces n'{q) to be a multiple 
of n{q), therefore zero, which makes the points of Cq planar and is a contradiction. □ 



4 Straight Creases Stay Straight 

Next we show that straight edges of a crease pattern must actually fold to straight line segments 
in space. 

Theorem 10 If ^ is a geodesic crease in a creased flat surface M with fold angle distinct from 
±180°, then -) is a segment in M^. 

Proof: The creased surface M decomposes by definition into a complex of uncreased surfaces, 
creases, and vertices. A point p in the interior of 7 is therefore on the boundary of two uncreased 
pieces; call them S and T. Let Sp and Tp be the tangent planes to S and T respectively at p. 
Because 7 is by hypothesis not a proper semicrease, has no semivertices along it, and has a fold 
angle distinct from ±180°, there is some p G 7 where Sp 7^ Tp. By continuity, the same is true for 
a neighborhood in 7 of p; let U be the maximal such neighborhood. 

Now parametrize 7 by arclength and let p = 7(5). At each q = ^{t), the tangent vector "f'{t) 
lies in the intersection Sq 7^ T^; in U , this determines ^'{t) up to sign. Because S and T are C^, the 
tangent planes Sq and Tq are C^, hence so is 7'(t), and the curvature 7"(t) exists and is continuous. 

Now around any q £ U project 7 onto the tangent plane Sq. Because 7 is a geodesic, we get 
a curve of zero curvature at q, so 7"(t) must be perpendicular to Sq. Similarly 7"(t) ± Tq. But 
certainly -f"{t) ± -f'{t). So -f"{t) = 0. 

We have 7"(t) = for t in a neighborhood of s, so 7 is a segment on U. Further, by the 
considerations of Theorem [H the tangent planes Sq and Tq are constant on U. Therefore they 
remain distinct at the endpoints of U, and because U is maximal, these must be the endpoints of 
7 and 7 is a segment. □ 

Combining the previous two theorems, we deduce that polygonal faces of the crease pattern 
with no boundary edges must indeed stay planar: 

Corollary 11 If an uncreased region of a creased flat surface M is piecewise geodesic and entirely 
interior to M , then the region lies in a plane. 
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5 Nonexistence of Pleated Hyperbolic Paraboloid 



Now we can apply our theory to prove nonfoldability of crease patterns. First we need to formally 
define what this means. 

Definition 4 A piece of paper is a planar compact 2-manifold. A crease pattern is a graph em- 
bedded into a piece of paper, with each edge embedded as a non-self-intersecting curve. A proper 
folding of a crease pattern is an isometric embedding of the piece of paper into 3D whose image is 
a good surface such that the union of vertices and edges of the crease pattern map onto the union 
of vertices and creases of the good surface. A rigid folding is a proper folding that maps each face 
of the crease pattern into a plane (and thus acts as a rigid motion on each face). 

Note that a proper folding must fold every edge of the crease pattern by an angle distinct from 
(to be a crease) and from ±180° (to be an embedding). We call such fold angles nontrivial. Also, 
one edge of a crease pattern may map to multiple creases in 3D, because of intervening semivertices. 

The key property we need from the theory developed in the previous sections is the following 
consequence of Corollary [TTl 

Corollary 12 For any crease pattern made up of just straight edges, any proper folding must fold 
the interior faces rigidly. 

We start by observing that the center of the standard crease pattern for a pleated hyperbolic 
paraboloid has no proper folding. 

Lemma 13 Any crease pattern containing four right triangular faces, connected in a cycle along 
their short edges, has no rigid folding. 

Proof: This well-known lemma follows from, e.g., |DQ02l Lemma 9]. For completeness, we give a 
proof. Let f i , v^, and denote the direction vectors of the four short edges of the triangular 
faces, in cyclic order. By the planarity of the faces, the angle between adjacent direction vectors is 
kept at 90°. Thus the folding angle of edge i equals the angle between Vi-i and ViJ^i (where indices 
are treated modulo 4) . If edge 2 is folded nontrivially, then vi and f 3 are nonparallel and define a 
single plane 11. Because V2 is perpendicular to both vi and ^3, V2 is perpendicular to 11. Similarly, 
Vi is perpendicular to 11. Thus V2 and f 4 are parallel, and hence edge 3 is folded trivially. Therefore 
two consecutive creases cannot both be folded nontrivially. □ 



Corollary 14 The standard crease pattern for a pleated hyperbolic paraboloid (shown in Fig- 
l(a)^ , with n > 2 rings, has no proper folding. 



ure 



Proof: With n > 2 rings, the four central triangular faces are completely interior. By Corollarv[T2l 
any proper folding keeps these faces planar. But Lemma [T3l forbids these faces from folding rigidly. 

□ 

The standard crease pattern for a pleated hyperbolic paraboloid cannot fold properly for a 
deeper reason than the central ring. To prove this, we consider the holey crease pattern in which 



the central ring of triangles has been cut out, as shown in Figure 3(a) If there were n rings in the 



initial crease pattern (counting the central four triangles as one ring), then n — 1 rings remain. 



Theorem 15 The holey crease pattern for a pleated hyperbolic paraboloid (shown in Figure 3(a)), 
with n — 1 > 3 rings, has no proper folding. 
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(a) Holey mountain-valley pattern for the 
pleated hyperbolic paraboloid. 



(b) Holey concentric pleat mountain-valley pattern. 



Figure 3: Holey mountain-valley patterns which have no proper foldings. 

Proof: Consider any nonboundary square ring of the crease pattern. By Corollary [T2l the four 
trapezoidal faces each remain planar. Any folding of these four faces in fact induces a folding of 
their extension to four meeting right triangles. But Lemma [13] forbids these faces from folding 
rigidly. □ 

A different argument proves nonfoldability of a more general pleated crease pattern. Define the 
concentric pleat crease pattern to consist of n uniformly scaled copies of a convex polygon P, in 
perspective from an interior point p, together with the "diagonals" connecting p to each vertex of 
each copy of P. The outermost copy of the polygon P is the boundary of the piece of paper, and in 
the holey concentric pleat mountain-valley pattern we additionally cut out a hole bounded by the 
innermost copy of P. Thus n — 1 rings remain; Figure |3(b)| shows an example. 

First we need to argue about which creases can be mountains and valleys. 

Definition 5 A mountain-valley pattern is a crease pattern together with an assignment of signs 
(+1 for "mountain" and —1 for "valley") to the edges of a crease pattern. A proper folding of a 
mountain-valley pattern, in addition to being a proper folding of the crease pattern, must have the 
signs of the fold angles (relative to some canonical orientation of the top side of the piece of paper) 
match the signs of the mountain-valley assignment. 

Lemma 16 A single-vertex mountain-valley pattern consisting of entirely mountains or entirely 
valleys has no proper rigid folding. 

Proof: If we intersect the piece of paper with a (small) unit sphere centered at the vertex, we 
obtain a spherical polygon whose edge lengths match the angles between consecutive edges of the 
crease pattern. (Here we rely on the fact that the vertex is intrinsically flat, so that the polygon lies 
in a hemisphere and thus no edges go the "wrong way" around the sphere.) The total perimeter of 
the spherical polygon is 360°. Any rigid folding induces a non-self-intersecting spherical polygon, 
with mountain folds mapping to convex angles and valley folds mapping to reflex angles, or vice 
versa (depending on the orientation of the piece of paper). To be entirely mountain or entirely 
valley, the folded spherical polygon must be locally convex, and by non-self-intersection, convex. 
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But any convex spherical polygon (that is not planar) has perimeter strictly less than 360° |Hal85l 
page 265, Theorem IV], a contradiction. □ 



Lemma 17 Consider a rigidly foldable degree-A single-vertex mountain-valley pattern with angles 
^1; ^2; ^3, o-nd 6i in cyclic order. Then exactly one edge of the mountain-valley pattern has sign 
different from the other three, and if 62 + 0^ > 180°, then the unique edge cannot be the one between 
62 and 63. 

Proof: Again we intersect the piece of paper with a (small) unit sphere centered at the vertex to 
obtain a spherical polygon, with edge lengths 9i, 62, O3, and 6*4, and whose convex angles correspond 
to mountain folds and whose reflex angles correspond to valley folds, or vice versa. By Lemma [TBI 
at least one vertex is reflex, and thus the remaining vertices must be convex. (The only non-self- 
intersecting spherical polygons with only two convex vertices lie along a line, and hence have no 
reflex vertices. Here we rely on the fact that the vertex is intrinsically flat, so that the polygon lies 
in a hemisphere, to define the interior.) The two edges incident to this vertex form a triangle, by 
adding a closing edge. The other two edges of the quadrilateral also form a triangle, with the same 
closing edge, that strictly contains the previous triangle. The latter triangle therefore has strictly 
larger perimeter than the former triangle, as any convex spherical polygon has larger perimeter 
than that of any convex spherical polygon it contains |Hal85|. page 264, Theorem III]. The two 
triangles share an edge which appears in both perimeter sums, so we obtain that the two edges 
incident to the reflex angle sum to less than half the total perimeter of the quadrilateral, which is 
360°. Therefore they cannot be the edges corresponding to angles 62 and ^3. □ 

Now we can prove the nonfoldability of a general concentric pleat: 

Theorem 18 The holey concentric pleat crease pattern (shown in Figure \3(b)^ , with n — 1 > 4 
rings, has no proper folding. 

Proof: First we focus on two consecutive nonboundary rings of the crease pattern, which by 
Corollary [12] fold rigidly. Each degree-4 vertex between the two rings has a consecutive pair of 
angles summing to more than 180° (the local exterior of P), and two consecutive pairs of angles 
summing to exactly 180° (because the diagonals are collinear). By Lemma [T71 the interior diagonal 
must be the unique crease with sign different from the other three. Thus all of the creases between 
the rings have the same sign, which is the same sign as all of the diagonal creases in the outer ring. 

Now focus on the outer ring, whose diagonal creases all have the same sign. Any folding of the 
faces of a ring in fact induces a folding of their extension to meeting triangles. (In the unfolded 
state, the central point is the center p of scaling.) Thus we obtain a rigid folding of a crease pattern 
with a single vertex p and one emanating edge per vertex of P, all with the same sign. But such a 
folding contradicts Lemma \W\ . □ 



6 Existence of Triangulated Hyperbolic Paraboloid 

In contrast to the classic hyperbolic paraboloid model, we show that triangulating each trapezoidal 
face and retriangulating the central ring permits folding: 

Theorem 19 The two mountain-valley patterns in Figure\^ with mountains and valleys matching 



the hyperbolic paraboloid of Figure 1(a) , have proper foldings, uniquely determined by the fold angle 
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9 of the central diagonal, that exist at least for n = 100 and 9 G {2°, 4°, 6°, . . . , 178°} for the 
alternating asymmetric triangulation, and at least for n and 9 shown in TableUlfor the asymmetric 
triangulation. For each 6 S {2°, 4°, . . . , 178°}, the asymmetric triangulation has no proper folding 
for n larger than the limit values shown in Table d 





(a) Asymmetric triangulation. 



(b) Alternating asymmetric triangulation. 



Figure 4: Two foldable triangulations of the hyperbolic paraboloid crease pattern (less one diagonal 
in the center). 



9 


2° 


4° 


6° 


8° 


10° 


12° 


14° 


16° 


18° 


20° 


22° 


24° 


26° 28° 


30° 
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176° 


178° 


n 


9 


7 


7 


7 


7 


7 


7 


5 


5 




5 


3 


3 


3 


3 



Table 1: The largest n for which the asymmetric triangulation has a proper folding, for each 
G {2°, 4°, . . . , 178°}. (By contrast, the alternating asymmetric triangulation has a proper folding 
for n = 100 for all such 9.) Interestingly, for 9 not too large, n • ^ is roughly 270°. 



Proof: The proof is by construction: we give a construction which implies uniqueness, and then 
use the resulting algorithm to construct the exphcit 3D geometry using interval arithmetic and a 
computer program. 

To get started, we are given the fold angle 9 between the two triangles of the central square. 
By fixing one of these triangles in a canonical position, we obtain the coordinates of the central 
square's vertices by a simple rotation. 

We claim that all other vertices are then determined by a sequence of intersection-of-three- 
spheres computations from the inside out. In the asymmetric triangulation of Figure 4(a), the 
lower-left and upper-right corners of each square have three known (creased) distances to three 
vertices from the previous square. Here we use Theorem [10] which guarantees that the creases 
remain straight and thus their endpoints have known distance. Thus we can compute these vertices 
as the intersections of three spheres with known centers and radii. Afterward, the lower-right and 
upper-left corners of the same square have three known (creased) distances to three known vertices. 
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one from the previous square and two from the current square. Thus we can compute these vertices 



also as the intersection of three spheres. In the alternating asymmetric triangulation of Figure 4(a) 



half of the squares behave the same, and the other half compute their corners in the opposite order 
(first lower-right and upper-left, then lower-left and upper-right). 

The intersection of three generic spheres is zero-dimensional, but in general may not exist and 
may not be unique. Specifically, the intersection of two distinct spheres is either a circle, a point, 
or nothing; the further intersection with another sphere whose center is not collinear with the first 
two spheres' centers is either two points, one point, or nothing. When there are two solutions, 
however, they are reflections of each other through the plane containing the three sphere centers. 
(The circle of intersection of the first two spheres is centered at a point in this plane, and thus the 
two intersection points are equidistant from the plane.) 

For the hyperbolic paraboloid, we can use the mountain- valley assignment (from Figure 1(a)) 
to uniquely determine which intersection to choose. In the first intersection of three spheres, one 
solution would make two square creases mountains (when the solution is below the plane) and 
the other solution would make those creases valleys. Thus we choose whichever is appropriate for 
the alternation. In the second intersection of three spheres, one solution would make a diagonal 
crease mountain, and the other solution would make that crease valley. Again we choose whichever 
is appropriate for the alternation. Therefore the folding is uniquely determined by 9 and by the 
mountain-valley assignment of the original hyperbolic paraboloid creases. 

This construction immediately suggests an algorithm to construct the proper folding. The 
coordinates of intersection of three spheres can be written as a radical expression in the center 
coordinates and radii (using addition, subtraction, multiplication, division, and square roots). See 
|Wik08j for one derivation; Mathematica's fully expanded solution for the general case (computed 
with Solve) uses over 150,000 leaf expressions (constant factors and variable occurrences). Thus, if 
the coordinates of the central square can be represented by radical expressions (e.g., is a multiple 
of 15°), then all coordinates in the proper folding can be so represented. Unfortunately, we found 
this purely algebraic approach to be computationally infeasible beyond the second square; the 
expressions immediately become too unwieldy to manipulate (barring some unknown simplification 
which Mathematica could not find). 

Therefore we opt to approximate the solution, while still guaranteeing that an exact solution 
exists, via interval arithmetic |Hay03 , [MKCOQl I AH83 j . The idea is to represent every coordinate x as 
an interval [xj;,,^;^] of possible values, and use conservative estimates in every arithmetic operation 
and square-root computation to guarantee that the answer is in the computed interval. For example, 
[a, b] + [c, d] = [a + c,b + d] and [a, b] • [c, d] = [min{a • c,a • d,b • c,b • d}, max{a ■ c,a ■ d,b ■ c,b ■ d}], 
while \/[a, b] requires a careful implementation of a square-root approximation algorithm such as 
Newton's Method. The key exception is that \/[a, b] is undefined when a < 0. A negative square 
root is the only way that the intersection of three spheres, and thus the folding, can fail to exist. If we 
succeed in computing an approximate folding using interval arithmetic without attempting to take 
the square root of a partially negative interval, then an exact folding must exist. Once constructed, 
we need only check that the folding does not intersect itself (i.e., forms an embedding). 

We have implemented this interval-arithmetic construction in Mathematica; refer to Appendix 1X1 
Using sufficiently many (between 1,024 and 2,048) digits of precision in the interval computations, 
the computation succeeds for the claimed ranges of n and 9 for both triangulations. Table [2] shows 
how the required precision grows with n (roughly linearly), for a few different instances. Figure [5] 
shows some of the computed structures (whose intervals are much smaller than the drawn line 
thickness). The folding construction produces an answer for the asymmetric triangulation even for 
n = 100 and 9 G {2°, 4°, . . . , 178°}, but the folding self-intersects for n larger than the limit values 
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shown in Table [TJ 



□ 



digits of precision 


16 


32 


64 


128 


256 


512 


1024 


2048 


nfoi 9 = r 


3 


6 


12 


22 


41 


76 


> 100 




n for 6* = 1° alt. 


3 


6 


12 


24 


43 


79 


> 100 




n for 6* = 45° alt. 


3 


5 


10 


18 


32 


58 


> 100 




n for (9 = 76° alt. 


2 


5 


9 


16 


29 


53 


95 


> 100 



Table 2: Number n of triangulated rings that can be successfully constructed using various preci- 
sions (measured in digits) of interval arithmetic. 



We conjecture that this theorem holds for all n and all 9 < 180° for the alternating asymmetric 
triangulation, but lack an appropriately general proof technique. If the construction indeed works 
for all 9 in some interval [0, 0), then we would also obtain a continuous folding motion. 




(a) Asymmetric triangulation, 6 = 8°, n = 16. (b) Alternating asymmetric triangu- 

lation, 6 = 30°, n = 16. 



Figure 5: Proper foldings of triangulated hyperbolic paraboloids. 



Interestingly, the diagonal cross-sections of these structures seem to approach parabolic in the 
limit. Figure [6] shows the x = y >0 cross-section of the example from Figure 5(b), extended out to 
n = 100. The parabolic fit we use for each parity class is the unique quadratic polynomial passing 
through the three points in the parity class farthest from the center. The resulting error near the 
center is significant, but still much smaller than the diagonal crease length, Least-square fits 
reduce this error but do not illustrate the limiting behavior. 



7 Smooth Hyperbolic Paraboloid 

Given a smooth plane curve F and an embedding of F in space as a smooth space curve 7, previous 
work |FT99j has studied the problem of folding a strip of paper so that a crease in the form of F in 
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40 



50 



(b) Absolute difference: fit minus actual. 




10 20 30 40 

(c) Relative difference: fit over actual. 



r 5 10 15 20 

(a) Actual zig-zag and parabolic fits. 

Figure 6: Planar cross-section of alternating asymmetric triangulation, 9 = 30°, n = 100, with 
parabolic fits of each parity class based on the last three vertices. 
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the plane follows the space curve 7 when folded. The main theorem from this work is that such a 
folding always exists, at least for a sufficiently narrow strip about T, under the condition that the 
curvature of 7 be everywhere strictly greater than that of T. 

Further, with some differential geometry described in |FT99j . it is possible to write down exactly 
how the strip folds in space; there are always exactly two possible choices, and additionally two 
ways to fold the strip so that T lies along 7 but remains uncreased. 

Based on some preliminary work using these techniques, we conjecture that the circular pleat 
indeed folds, and that so too does any similar crease pattern consisting of a concentric series of 
convex smooth curves. Unfortunately a proof remains elusive. Such a proof would be the first 
proof to our knowledge of the existence of any curved-crease origami model, beyond the local 
neighborhood of a single crease. 
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A Mathematica Construction of Triangulated Hyperbolic Paraboloid 

SW = 1; SE = 2; NE = 3; NW = 4; 
next [d_] := 1 + Mod[d, 4] ; 
prev[d_] := 1 + Mod[d + 2, 4]; 

square [i_] := {{-i, -i}, {i, -i}, {i, i}, {-i, i}}; 
side[i_] := 2 i; 

diagonal [i_] := Sqrt [2] side[i]; 
diagpiece [i_] := Sqrt [2] ; 

cross[i_] := Sqrt[(i - (i - 1))~2 + (i + (i - 1))~2]; 
n = 100; 

bigside := side[n]; 
bigdiagonal : = diagonal [n] ; 
trapezoid [i_ , d_] := 

square [i] [ [{d, nextCd]}]] "Join" squared - 1] [ [{next [d] , d}] ] ; 

(* Note: \[Theta\] is half the fold angle. *) 
algebraicmiddle := { 

{-1, -1, 0>, {-Cos [\ [Theta] ] , Cos [\ [Theta] ] , Sqrt [2] Sin [\ [Theta] ] > , 

{1, 1, 0}, {Cos [\ [Theta] ] , -Cos [\ [Theta] ] , Sqrt [2] Sin [\ [Theta] ]» ; 
precision = 1024; 

middle := Map [Interval [N [#, precision]] &, algebraicmiddle, {2}]; 

orientation [a_ , b_, c_, d_] := 

Sign [Det [Append [#, 1] & /@ {a, b, c, d}] ] ; 
pickorientation [choices_ , a_, b_, c_, pos_] := 
Module [ 

{select = Select [choices , 

If [pos, orientation [a, b, c, #] >= &, 

orientation [a, b, c, #] <= &]]}, 

If [Length [select] > 1, 

Print ["Warning: Multiple choices in pickorientation..."]]; 
If [Length [select] > 0, select [[1]], 
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Print ["Couldn't find mountain/valley choice"]]]; 
triangles [1] := {triangulated [ [1 , {1, 2, 3}]], triangulated [ [1 , {3, 4, 1}]]}; 

threespheres = 

Solve [{SquaredEuclideanDistance [{xl , yl, zl}, {x, y, z}] == dl~2, 
SquaredEuclideanDistance [{x2 , y2, z2}, {x, y, z}] == d2~2, 
SquaredEuclideanDistance [{x3 , y3, z3}, {x, y, z}] == d3~2}, {x, y, z}] ; 
solvethreespheres[{Xl_, Yl_, Zl_}, Dl_, {X2_, Y2_, Z2_}, 
D2_, {X3_, Y3_, Z3_}, D3_] := 
threespheres /. {xl -> XI, yl -> Yl, zl -> Zl, x2 -> X2, y2 -> Y2, 

z2 -> Z2, x3 -> X3, y3 -> Y3, z3 -> Z3, dl -> Dl, d2 -> D2, d3 -> D3}; 

checklengths := Module [{i, j, t, tri, error, maxerror, who}, 
maxerror = ; 

who = "I did not see any error!"; 
For[i =1, i <= Length [triangulated] , i++. 
For [t = 1 , t <= Length [triangles [i] ] , t++, 
tri = triangles [i] [ [t] ] ; 
For[j = 1, j <= 3, 
error = 

Min[Abs [EuclideanDistance [tri [[j]] , tri [ [Mod [j, 3] +1]]] - #] & /O 
If [i == 1, {side[i], 2 diagpiece [i] } , 

{side[i], side[i - 1], cross [i], diagpiece [i] }]] ; 
If [error >= maxerror, 

who = {"i=", i, " triangle ", t, " edge ", j, "-", 

Mod[j, 3] +1, " has error ", error, " = ", N [error], ": ", 
EuclideanDistance [tri [ [j] ] , tri [ [Mod [j, 3] + 1]]], 
" vs. ", {side[i], side[i - 1], cross [i], diagpiece [i] }} ; 
maxerror = error]]]]; 
Print m who] ; 

(* Collision detection *) 

stL = Solve [{(1 - L) ax + L bx == 

px + s (qx - px) + t (rx - px) , (1 - L) ay + L by == 
py + s (qy - py) + t (ry - py) , (1 - L) az + L bz == 
pz + s (qz - pz) + t (rz - pz)}, {s, t, L}] [[1]] ; 

pierces [seg_ , tri_] := Module [{parms} , 
Check [Quiet [ 

parms = stL / . {ax -> seg [ [1] ] [ [1] ] , ay -> seg [ [1] ] [ [2] ] , 

az -> seg[[l]] [[3]] , bx -> seg [ [2] ] [ [1] ] , by -> seg [ [2] ] [ [2] ] , 

bz -> seg[[2]] [[3]] , px -> tri[[l]] [[1]] , 
py -> tri[[l]] [[2]] , pz -> tri[[l]] [[3]] , qx -> tri [ [2] ] [ [1] ] , 

qy -> tri[[2]] [[2]] , qz -> tri [ [2] ] [ [3] ] , 
rx -> tri[[3]] [[1]] , ry -> tri [ [3] ] [ [2] ] , 
rz -> tri[[3]] [[3]]}; 
spt = (s + t) /. parms; LL = L / . parms; 
Return [(s >0&&t>0&&s+t<l&&0<L<l) /. 
parms], {Power :: "infy" , \ [Infinity] :: "indet"}] , 
Return [False] , {Power :: "infy" , \ [Infinity] :: "indet"}] ] ; 
checkcollision := Module [{i, j, t, t2, tri, seg}, 
For[i =1, i <= Length [triangulated] - 1, i++, 
PrintTemporary ["Ring ", i, "vs. ", i + 1] ; 
For [t = 1 , t <= Length [triangles [i] ] , t++. 
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tri = triangles [i] [ [t] ] ; 

For[t2 =1, t2 <= Length [triangles [i +1]], t2++, 
For[j = 1, j <= 3, 
seg = triangles [i + 1] [ [t2] ] [ [{j , Mod[j, 3] + 1>]]; 
If [Length[Intersection[seg, tri]] > 0, Continue[]]; 
If [pierces [seg, tri], 

Print ["COLLISION! ", i, t, " vs. ", i + 1, t2, j] ; 

Return [{seg , tr i}] ] ] ] ] ] ; 
Print["No collision."]]; 

adaptiveprecisiontest [mint_ , maxt_, tstep_: 1, minprecision_ : 8, maxprecision_ : 65536] := 
Module [{t}, 

precision = minprecision; 
Print ["n: " , n] ; 

For [t = mint, t <= maxt , t = t + tstep, 
Print [t, " degrees:"]; 
\[Theta] = t \ [Pi] /ISO; 

For [True, precision <= maxprecision, precision = 2 precision. 
Print ["precision: ", precision]; 
If [computetriangulated == n. 

Print ["Required precision for ", t, " degrees:", precision]; 

checklengths ; checkcollision; Break[]]]]]; 

Print ["ALTERNATING ASYMMETRIC TRIANGULATED FOLDING"]; 

oddoutset [i_, d_ , last_] := 

pickorientation[{x, y, z} / . # & /@ 

solvethreespheres [last [ [d] ] , diagpiece [i] , last [ [next [d] ]] , 

cross [i], last [ [prev [d] ] ] , cross [i]], 
last[[prev[d]]] , last[[d]], last [ [next [d] ]] , Oddq[i]]; 
evenoutset [i_ , d_ , last_, odds_] : = 
pickorientation[{x, y, z} / . # & /@ 

solvethreespheres [last [ [d] ] , diagpiece [i] , odds [ [next [d] ] ] , 

side [i] , odds [ [prev [d] ] ] , side [i] ] , 
odds [[prev [d]]] , last[[d]], odds [ [next [d] ]] , OddQ[i]]; 
computetriangulated := Module [{i, odds}, 
triangulated = {middle}; 
For[i =2, i <= n, i++, 
If [OddQ[i] , 

last = RotateRight [triangulated [ [-1] ]] , 
last = triangulated [ [-1] ]] ; 
odds = Table [If [OddQ [d] , oddoutset [i, d, last]], {d, 1, 4}]; 
Print ["Round ", i, " odds done"]; 
If [Count [odds. Null] > 2, Break []]; 
both = 

Table [If [OddQ [d] , odds[[d]], evenoutset [i , d, last, odds]], {d, 1, 4}]; 
If [OddQ [i] , 

AppendTo [triangulated, RotateLef t [both] ] , 

AppendTo [triangulated, both]]; 
Print ["Round ", i, " evens done"]; 
If [Count [triangulated [[-1]] , Null] > 0, Break []]]; 
Print [i - 1, " rounds complete"]; 

i-1]; 



18 



triangles [i_] := 
Flatten [Table [ 

If[OddQ[d] == EvenQ[i], 

{{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [ [i - 1, prev[d]]]}, 
{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [ [i - 1, next[d]]]}}, 
{{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [ [i , prev[d]]]}, 
{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [[i, next[d]]]»], {d, 1, 4>] , 1]; 

adaptiveprecisiontest [1 , 90]; 

Print ["ASYMMETRIC TRIANGULATED FOLDING"] ; 

computetriangulated := Module [{i, odds}, 
triangulated = {middle}; 
For[i = 2, i <= n, i++, 

odds = Table[If [OddQ[d] , oddoutset[i, d, triangulated [ [-1] ]]] , {d, 1, 4}]; 
Print ["Round ", i, " odds done!"]; 
If [Count [odds , Null] > 2, Break []]; 
AppendTo [triangulated, 

Table [If [OddQ [d] , odds [ [d] ] , 

evenoutset [i , d, triangulated [ [-1] ] , odds]], {d, 1, 4}]]; 
Print ["Round ", i, " evens done"]; 
If [Count [triangulated [[-!]] , Null] > 0, Break []]]; 
Print [i - 1, " rounds complete"]; 

i-1]; 

triangles [i_] := 
Flatten [Table [ 
If [OddQ [d] , 

{{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [ [i - 1, prev[d]]]}, 
{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [ [i - 1, next[d]]]}}, 
{{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [ [i , prev[d]]]}, 
{triangulated [ [i , d] ] , triangulated [ [i - 1, d] ] , 
triangulated [[i, next [d] ] ] }}] , {d, 1, 4}], 1] 

adaptiveprecisiontest [1 , 90]; 
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